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Abstract 

Context. A small number of young stellar objects show signs of a halo-like structure of optically thin dust, in addition to a circum- 
stellar disk. This halo or torus is located within a few AU of the star, but its origin has not yet been understood. 
Aims. A dynamically excited cloud of planetesimals colliding to eventually form dust could produce such a structure. The cause of 
the dynamical excitation could be one or more planets, perhaps on eccentric orbits, or a migrating planet. This work investigates an 
inwardly migrating planet that is dynamically scattering planetesimals as a possible cause for the observed structures. If this mecha- 
nism is responsible, the observed halo-like structure could be used to infer the existence of planets in these systems. 
Methods. We present analytical estimates on the maximum inclination reached owing to dynamical interactions between planetesi- 
mals and a migrating planet. In addition, a symplectic integrator is used to simulate the elfect of a migrating planet on a population 
of planetesimals. Collision time scales are estimated for the resulting population of planetesimals and the size distribution of the dust 
created in catastrophic collisions is determined. 

Results. It is found that an inwardly migrating planet is only able to scatter the material it encounters to highly-inclined orbits if the 
material is on an eccentric orbit. Such eccentric orbits can be the result of resonance trapping and eccentricity pumping. Simulations 
show that for a certain range of migration rates and planet masses, resonance capture combined with planetary migration indeed 
causes the planetesimals to reach eccentric orbits and subsequently get scattered to highly-inclined orbits. The size distribution of the 
resulting dust is calculated determined to find the total mass and optical depth, which are found to compare reasonably well with the 
observed structures. 

Conclusions. Dynamical scattering of planetesimals caused by a planet migrating in, followed by the grinding down of these plan- 
etesimals to dust grains, appears to be a promising explanation for the inferred circumstellar dust clouds. Further study is needed to 
see if the haloes can be used to infer the presence of planets. 



Key words. Protoplanetary disks - Planet-disk interactions - Planets and satellites: dynamical evolution and stability - Stars: circum- 
stellar matter 
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1 . Introduction and on the temperatures and particles found in th e dilferent parts 

ly-j of the circumstellar disk (iDullemond et al.l[2007h . 

O The vast majority of stars are bom with a circumstellar disk. Over the past two decades, studying SEDs has improved 

T-H Observations show that in young clusters the percentage of stars our knowledge of circumstellar disks significantly and has made 

^ , with disks is close to 10 0%. The typical life times for disks ap- us aware of structur es as flaring disks, puffed up inner rims, 

■ pear to be ai'ound 6 Mvr (iHaisch e t al. 2001 ). Circumstellar disks gaps, and holes (e.g. IDullemond et all 1200 ll l2007h . Holes and 

• ^ around young stars are believed to be the sites of planet forma- gaps are conside red to be possible signposts of planet formation 

tion. Solids in the disk grow to form larger objects, planetes- (I Armitag3l20 1 Oh . 

H imals, and eventually planets. Eventually, after the gas disap- Recent studies of (pre-)transitional disks indicate a possible 

. . . pears, a planetary system will remain, containing planets and, if new structure. Several young stellar objects (YSOs) appear to 

there is room in the dynamical structure, a debris belt. show signs of an optically thin population of dust with a large 

„, . J. , 1 , J- 1 , . • . scale height close to the star These structures appear to be opti- 

Observations of protoplanetary disks at ever higher spa- n i • i i n i i j ■ r 

. , , . , ^ , ■' . . ... . J. cally thin, but they still absorb and reprocess a significant frac- 

tial resolution have provided promising indirect indications ■-•,,'„,•,,,•, • .i j . i j- -i j 

„ , \. . ^ , — ; — — — —n rrrm tion of the stellar light, which requires the dust to be distributed 

of ongoing planet formation (e.g. IThalma nn et al.l I2010t . ^ , u- u i u • u. j . j . 

r; 7T — p — r — nr^Trrri » 1 • j- \- c m a structure with high scale height, and not compressed to- 

IJang-Condell & Kuchnen 120101) . Another strong indication of , ^, . , i rfL j- i o u i u • u. r.u . 

~; — ^ z — ', 1 — \: , J • 1 ward the midplane of the disk. Such a large height of the struc- 

the presence of planets comes from structure observed in the . . .• , , ■ ■ j • i-i ■ 

J Al . ■ 11 ■ i-i 1 • ■ J 1 • J- 1 ,,,r i — r~n tures IS puzzling, since disks that are in hydrostatic equilibnum 

dust that IS still visible a t later times in debns disks (Wiln er et al.l . ui. ^ . i . .u . n? — o tt ^ — IFTTT^ 

However, since most of the disks are difficult to resolve, stud- ity that the large height of the halo-like structure may be reached 

ies also focus on understanding the system's spectral energy dis- in a dynamical way by scattering planetesimals. This idea is in- 

tribution (SED). Comparing the shape of the SED to radiative spired by the fact that planetesimals in the Solar System (in the 

transfer models gives information on the geometry of the disk Asteroid Belt and the Kuiper Belt) can have substantial incli- 
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nations, caused by interactions with planets. If this mechanism 
is found to be capable of explaining the observed optically thin 
structures, their presence could be used to infer the existence of 
planets in these systems. This work will focus on the question 
of whether a inwardly migrating planet could be responsible for 
creating such an optically thin halo-like structure. 

Section|2]discusses three YSOs that show this new structure. 
Section |3] discusses high inclinations in our own Solar System, 
and quantitatively formulates a mechanism that could cause 
highly-inclined orbits in a system with one planet. In section 
|4] numerical simulations of the effect of a migrating planet on 
a planetesimal population are presented. The grinding down of 
the resulting planetesimal population to form small dust grains 
is covered in section|5] Finally, a discussion and conclusions are 
formulated in sections |6] and |7] 

2. Optically Thin Dust 

Recently, a number of transitional disks have indications of pos- 
sessing optically thin dust with large scale height: LkCa 15, HD 
142527 and HD 163296. These systems are discussed below. 

2.1. LkCa 15 

LkCa 15 is a low-mass pre-main-sequenc e star located at abou t 
140 pc in the Taurus star-forming region (lEspaillat et al.ll2007l) . 
lEspaillat et al.l (l2007l) compared different models to SEDs ob- 
tained with the Infra-red Array Camera (IRAC) aboard Spitzer 
to study LkCa 15. An optically thick outer disk is found to have 
an inner radius of 46 AU, containing about O.lMp. The i nner 
region of the disk is a bit more complex and Espail lat et al.l pro- 
vide two options which could explain the shape of the SED in 
the near-IR: an optically thick disk from 0.12 to 0.15 AU com- 
bined with 4 X 1O""M0 of optically thin dust between 0.15 and 
4 AU, or, 5 X 1 0" " Mq of optically thin dust between 0.12 and 4 
AU. The first option is said to give the better fit. Both scenarios 
include an optically thin part, but it is stressed no dust can exist 
further out tha n ~ 5 AU sin ce the contribution at 20 fim becomes 
too strong (Espaillat et al.ll2007ll) . 

In 2009, Mulders et al.l investigate the inner region of LkCa 
15 in more detail and find an optically thick inner disk will in- 
fluence the outer disk as well. Because the optically thick inner 
disk will obscure the outer disk, the outer disk needs to be blown 
up in the vertical direction in order for it to fit the observed SED. 
Alternatively, the inner disk could be optically thin, but this is 
only possible if the material i s distributed over a la rge solid an- 
gle, virtually making it a halo. Mulders et al. (2010) used the 2D 
radiative transfer code MCMax (Min et al. 2009) to model the 
observed SED. Mulde rs et al.l conclude it is hard to explain the 
large scale height of the inner dust material. One offered expla- 
nation is the violent scattering and disruption of planetesimals, 
which would end up in them colhding and filling the inner region 
wit h dust. 

lEspaillat et al.l ( l2OI0l) use new observations of the near-IR 
excess and find it to have the spectral shape of a blackbody, lead- 
ing them to prefer the option of an optically thick inner disk and 
a blown-up outer disk. The inner disk is found to lie within 0. 19 
AU and the upper limit for the mass it contains is found to be 
2 X 10'* Mq. 

Whether the inner disk of LkCa 15 contains an optically thin 
dust component remains a matter of debate, but if optically thin 
dust is responsible for the near-IR flux, the scale height of the 
material has to be very large. 



2.2. HD 142527 

Another interesting system is t hat of HD 14 2527. The star was 
categorized as an F6III star b y Houk (1978) and later classified 
as a Herbig star (iWaelkens et al. 1996) . HD 142527 is about 
20 times as luminous as the sun, has a mass of 2.2 ± O.SM©, 
and shows a huge IR-flux (Fir - 0.92^*) which indicates 
t here must be a cons iderable amount of circumstellar material 
(IVerhoeff^et al.ll20Ilh . 

IVerhoeff^ et al.l used MCMax to model the circumstellar disk 
and find a model which accurately describes Spitzer and ISO IR 
spectra. A complicated structure is found for the system; an in- 
ner and an outer disk as well as an optically thin halo-like struc- 
ture. The inner disk stretches from 0.3 to 30 AU and has a puffed 
up inner rim which shadows the rest of the disk. The small dust 
grains in the inner disk have a total mass of 2.5 x 10 ''Mq. It is ar- 
gued most of the other material is already locked up in planetesi- 
mals or even planets. The degr ee of sedimentation a lso indicates 
an advanced state of evolution dVerhoefi^ et al.ll20TTl) . 

Then there is the optically thin halo, which stretches from 0.3 
to 30 AU and has an estimated dust mass of 1.3 x IO^'^Mq. The 
origin of this halo is unclear and again the dynamical scattering 
of planetesimals is offered as a speculative explanation. 

There is a gap in the disk stretching from 30 to 130 AU. 
IVerhoeff et al.l explain why photo-evaporation is unlikely to be 
the cause of this gap and speculates about the presence of Jupiter 
or Earth-like planets as the origin of the gap. Lastly there is the 
outer disk, located between 120 and 200 AU, and with a rim 
height of 60 AU. Such a large rim height can be reached in hy- 
drostatic equilibrium, but only if the disk is illuminated by the 
nearly unattenuated stellar radiation. A physical explanation for 
the high rim therefore requkes the warm dust to remain optic ally 
thin and, therefore, to be spread out to large scale heights. The 
total mass in small dust grains is estimated to be 1.0 x 10'^ Mq. 

2.3. HD 163296 

HD 163296 is an isolated Herbig Ae star located at about 122 
pc, which shows a NIR excess indicating a circumstellar disk 
Hillenbr and et al.l (fl992). In 2010, Benistv et al. present long- 
baseline spectro-interferometric observations using the AMBER 
instrument on VLTI. They find that the inner rim of a disk alone 
cannot reproduce their observations in the NIR and need an ad- 
ditional emission component between 0.10 and 0.45 AU, with a 
temperature of 1600 K and an optical depth of ~ 0.2. 

iBenistv et alj discusses whether hot gas could be held re- 
sponsible for the emission between 0.10 and 0.45 AU but con- 
clude that this is probably not the case. They argue that more 
detailed, non-LTE models that include the transition from opti- 
cally thin to optically thick layers in a dust-free environment are 
needed to completely rule out the possibility of hot gas being the 
main contributor to the NIR flux. 

The existence of refractory dust grains is explored as another 
possibility for the emission. An optically thin layer of dust grains 
is assumed to lie between an inner and outer radius, which are 
both inside the inner rim of the optically thick disk. The vertical 
optical depth of the dust is ass umed to be pr oportional to 1 /r, 
with r the distance to the star. iBenistv et al.l find the observed 
SED is reproduced best with the dust between 0. 1 and 0.45 AU. 
The total mass in small dust between is calculated and found to 
be 8.7 X 10"*^ Me for Graphite grains of 0.05 - 2jum, and 9.5 x 
10"^ Me for Iron grains of 0.02 - 2 //m in size. 
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3. Direct Scattering to Reach l-lighly-lnclined Orbits 

3.1. Highly-Inclined Orbits in our Solar System 

How do we build a population of planetesimals on highly- 
inclined orbits from a nearly flat transitional disk? Hints for a 
mechanism can be found in the Kuiper Belt in our own Solar 
System. The Kuiper Belt c an be split up into three dynamical 
classes (! Jewitt & Luul2000l) : The classical Kuiper Belt, the reso- 
nant Kuiper Belt, and the scattered Kuiper Belt. Classical Kuiper 
Belt Objects (KBOs) have low-eccentricity, low-inclination or- 
bits beyond 40 AU. The resonant KBOs are trapped in a mean 
motion resonance (MMR) with Neptune, and the scattered ones 
have probably suffered a close encounter with Neptune in the 
past and show a wide spread in inclination. Scattered KBOs tend 
to have eccentric orbits with large semi-major axes. Studies on 
the di stribution of inclinations of bodi es in these three popula- 
tions (lBrownll200"l1; iGulbis etal]|2010l) show they truly are dis- 
tinct classes and it is argued that the inclinations in the three 
classes have different origins. 

Observatio nal studies show resonant KBOs to have inclina- 
tions up to 30° (lBrownll200lHGulbis et al.ll20T^ . They are com- 
monly believed to have been trapped in an MMR with Neptune 
when the planet was migrating outward, having their eccentric- 
ity and inclination raised during the process. This mechanism 
explains Pluto's peculiar orbit (Malhotra 1993), and inclinations 
of up to ~ 15° in general for bodies outside Neptune's orbit 
(lMalhotralll995h . However, only a small fraction of the reso- 
nant KBOs in the sim ulations have t heir inclinations raised in 
excess of 10% (Malhotra etal.ll2000l) . Other mechanisms have 
been proposed to account for the o bserved inclina tions, includ- 
ing sweeping secular resonances (Li et al.' '200 81). and stirring 
by large planetesimals (Morbidelli & Valsecchil [l997l) . but no 
definitive answer has been given. 

Scattered KBOs show a wider spread in semi-major axes 
and their incl inations are usual ly greater than 10° with a peak 
around 20° (Gulb is et al.l l2010). It is this broad distribution in 
orbital parameters that indicates that their orbit s are probabl y 
the result of a close en counter with Neptune (iBrownl l2001b . 
iHahn & Malhotal d2005h conclude from simulations that about 
90% of the KBOs in the scattered disk might not have suffered a 
close encounter, but rather obtained their orbits from Neptune's 
resonances during the migration epoch. Their model is unable to 
reproduce the observed abundance of KBOs with / > 15° how- 
ever. 

From studying the Kuiper Belt, it appears there are 2 clear 
mechanisms that are able to increase an object's inclination 
drastically; resonance trapping and close encounters. The de- 
gree of inclination-raising by resonance trapping is determined, 
and therefore limited by, the migration distance of the planet. 
[Malh otra ( 1995) shows a migration Neptune over a distance of 
7 AU to cause an increase in eccentricity of 0.25 and in inclina- 
tion of about 10°. The inclination that can be reached by direct 
scattering is limited by the relative velocity between the planet 
and the object and the escape velocity at the planet's surface, as 
we shall see in section [32] It appears direct scattering is a more 
promising way to get bodies onto highly-inclined orbits, and it 
is this path that we will pursue in the rest of this work. 



on an orbit with negligible inclination is approaching the planet 
with an impact parameter b and a relative velocity v,ei. 

If b is small enough, the close encounter can be treated as a 
2 body interaction, without having to consider the central star. 
During this interaction, the magnitude of the relative velocity is 
cons erved but the direction is changed by a scattering angle of 
(e.g. IWeidenschillinglll975h 



6-2 arctan 



(1) 



We normalize distances to the radius of the planet and velocities 
to the Keplerian velocity 



V - — b - — 
vk ' ^pl ' 

This allows us to rewrite eq.[T]as follows 



0-2 arctan 



1 

2 \ Vi-el 



(2) 



(3) 



where Vesc = (2GMpi/7?pi)'^^ is the escape velocity at the surface 
of the planet. 

This equation shows that a scattering angle of 90° is reached 

for 



Vi-el 



1 

2b 



(4) 



Thus, for velocities Vrei > Vesc/ V2, an impact parameter smaller 
than the physical size of the planet is required. Such a planetesi- 
mal will crash into the planet. Due to gravitational focussing any 
planetesimal with an impact parameter smaller than 



1 + 



Vrel 



(5) 



will b e lost in a physical collision with the planet dSafronovl 
[19661). 

Thus, for a planet to be able to significantly scatter a plan- 
etesimal out of the orbital plane, the relative velocity of the en- 
counter needs to be about a factor 2 lower than the escape veloc- 
ity at the planet's surface. 

Suppose the planet is able to scatter the planetesimal out of 
the orbital plane by an angle of 90°. In the frame of the planet, 
the planetesimal's total velocity is now pointing perpendicular to 
the orbital plane. In the stellar frame however, the planetesimal 
is also moving with the planet at the Keplerian velocity, and we 
can calculate the new inclination using 



; - arctan 



Vrel 
V'K 



(6) 



This simple calculation shows that an inclination of 45° can be 
reached for Vrei - v'k. 

Summing up, planetesimals can only be scattered into 
highly-inclined orbits for a particular combination of Vrei, Vesc 
and vk. This condition can be written as 



3.2. Direct Planetesimal Scattering 

Suppose we have a planet of mass Mpi and radius Rpi orbiting 
a star of mass M,, on a circular orbit with Keplerian velocity 
Vk = -^jGMTJopi. We assume a planetesimal of mass m <k Mp\ 



Vk S Vrel ^ 



(7) 



We can write the relative velocity of the collision as a func- 
tion of orbital parameters of both objects. Recall that the planet 
moves on a circular orbit with semi-major axis Opi. We define the 
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orbit of the planetesimal by its eccentricity e, semi-major axis a, 
and inclination / relative to the planet. The relative velocity is 
then given by ( Weidenschilling 1975) 



= 3-2 / — (l-e2)cos/-— . 

flpi a 



(8) 



Recalling / = and taking a - flpi, we need an eccentricity of 
e > (3/4)'^^ » 0.87 to get Vrei > vk. Particles on orbits larger 
than that of the planet reach similar relative velocities for lower 
inclinations. However, if the particle is located inside the plane- 
tary orbit, for instance in the 1:2 MMR, it has a/opi ~ 0.63 and 
only reaches v'rei - vk for an eccentricity of 0.97. 

Since this study focusses on inward planetary migration, 
most of the material will have a < a^i. It is obvious from eq. 
[8] that this means that the planetesimals need to be on (very) ec- 
centric orbits. Section |43] will explore an explanation of these 
highly eccentric orbits. 

The remainder of this section will calculate post-encounter 
inclinations in more detail and study the conditions for which 
planetesimals are scattered out of the system on hyperbolic or- 
bits. For simplicity, we assume the orbits of the planetesimals to 
have e - \, causing the velocities of the planet and the planetes- 
imal to be perpendicular at the point of close encounter. Figure 
[1] shows the different quantities involved. Figure [TJ depicts the 
stellar frame. In this case the comet velocity points in the nega- 
tive x-direction, the planet moves in the direction of the positive 
y-axis, and the z-axis is perpendicular to the orbital plane. 

Figure [T]; shows the scattering in the frame of the planet. 
In this frame the planetesimal's velocity points in the negative 
x' -direction, and has a size Vrei determined by 



Vrel 



1 +vi 



(9) 



The shaded area depicts the plane in which the scattering over an 
angle of 6 occurs, which is the (y '-z)-plane. Unit vectors x and x' 
(and y and y ') lie in the orbital plane, but are rotated with respect 
to each other by an angle of 



a - arctan(vc). 



(10) 



as depicted in figure [TJ). 

As a planetesimal comes in (along the x '-axis) with impact 
parameter b, the scattering angle 9 is determined by eq. |3] We re- 
peat that a comet approaching with b < will have a physical 
collision with the planet. 

For the post-encounter comeQ inclination (/) and velocity 
(^post), it is not just the size of the impact parameter that matters, 
but also the direction of this offset with respect to the motion of 
the planet. To quantify this direction we use (p. If the planetesi- 
mal is aiming at a point slightly behind the planet {(p - n), it will 
be scattered in the planet's direction of motion, and therefore ac- 
celerated. If the planetesimal is aiming at a point directly above 
or below the planet {<p = ±7r/2), it is scattered out of the orbital 
plane. 

If we split the post-encounter velocity of the comet into a 
component in the z-direction and one in the orbital plane, Vpost.z 
post.xy respectively, we can calculate the total velocity to be 



(a) 



r 



'post 



(b) 




(c) 




Figure 1. (a) Scattering of a planetesimal in the stellar frame. 
The planetesimal comes in along the x-axis and the planet moves 
in the positive >'-direction. The z-axis is perpendicular to the or- 
bital plane, (b) Unit vectors x and x ' (and y and y ') lie in the 
orbital plane, but are rotated with respect to each other by an 
angle a. (c) The scattering in the frame of the planet. The plan- 
etesimal comes in along the x '-axis, with impact parameter b. It 
is scattered in the (y '-z)-plane over an angle 6. 



and the post-encounter incUnation using 



I - arctan 



^ post.z 



^ post.xy 



(12) 



Since the escape velocity from the central star is equal to V2vk, 
comets with Vpost > V2 will leave the system on a hyperbolic 
orbit. 

For a given combination of Vc, Vesc, b and (p, the resulting 
inclination and total velocity are uniquely defined and can be 
calculated using the following prescription: 

1. calculate Viei and bmin from eq.|9]and|5] 

2. determine angles a and 6, 

3. in the planet frame, rotate v,ei using 9 and (p, 

4. switch to the stellar frame by adding the planet's velocity 
under the correct angle a, 

to describe bodies of mass m « Mp,. Come/ is mostly used to refer to a 5. calculate Vpost using eq.[TT]and check if the particle is ejected, 
particle on an orbit with a large eccentricity. 6. find the comet's post-encounter inclination from eq. [12] 



and V 



^^post 



post,xy post,z' 



(11) 



The words comet and planetesimal are used throughout this paper 
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2 4 6 8 10 2 4 6 8 10 

h b 

Figure 2. The inclination of a comet after an encounter at Vrei = 1 .41 with a planet with f'esc, as a function of impact parameter and 
(f>. The white regions mark the combinations of b and where Vpost > V2. Regions where b < fomin are also coloured white. 



Figure |2] shows the post-encounter inclination of a comet 
coming in with Vc = 1 .41 as a function of b and for four differ- 
ent values of Vesc- The regions where the comet is lost because 
of ejection from the system or a collision with the planet are 
coloured white. 

For a given system, ejection is more likely to happen for 
comets with relatively small impact parameters, and (p ~ n. The 
latter is caused by these comets coming in behind the planet, 
and thus being accelerated in the stellar frame during the close 



encounter. The value of (p that is ideal for reaching high inclina- 
tions appears to lie between and 7r/2, because the comet not 
only gains velocity in the z-direction, but it is also slowed down 
in the xy-plane. 

The planet with the lowest mass can be seen not to be able 
to eject bodies, nor can it scatter them to orbits with inclinations 
> 5°. As the escape velocity at the planet's surface is increased, 
three things happen; 1) bmm increases, 2) the region where par- 
ticles are ejected grows larger (both in and in b), and 3) the 
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region where high inclinations are reached shifts towards higher 
impact parameters. Figure |2] confirms equation |7] as an escape 
velocity higher than the relative velocity is needed to scatter bod- 
ies to highly-inclined orbits. 

3.3. Eccentricity Pumping and Migration 

The relative velocity used in fig. |2]can only be reached for very 
eccentric planetesimal orbits. A possible explanatio n for high 
eccentricities is resonance trapping, see for instance iMalhotral 
([1993, 1995); Malhotra et al. (2000). In our own Solar System, 
the trapping and consequent dragging along of Plutinos in 
Neptune's 3:2 resonance, has raised their inclinations from 
nearly to up to 30° (see section im i. 

A planetesimal orbiting a star is in an MMR with a planet 
if their mean motions are locked , and t he resonant argument is 
librating (see iMurrav & Dermottl(ll999h for a detailed treatment 
of the resonant argument). The locations of the mean motion 
resonances of a planet orbiting a star at Upi are 



p + q 



2/3 



(13) 



with p and q integers, and the ratio of an MMR depicting the 
ratio of the orbital periods, so that a 1:2 resonance lies inside 
the planet's orbit while a 2:1 resonance lies outside. Resonances 
with \q\ = 1 are called first-order resonances and are often de- 
scribed with j instead of p. 

When the planet is migrating, the locations of the MMRs will 
change and bodies trapped within a resonance can be dragged 
along. This will not only change the dragged-along body's semi- 
major axis, but also various other orbital parameters. Malhotra 
(Il993h showed Pluto's peculiar orbit could be explained by it 
being stuck in Neptune's 3:2 resonance, as Neptune migrated 
out by ~ 7 AU. 

IWvattI (l2003h provides us with an expression for the ec- 
centricity of a planetesimal orbit, that has been dragged in a 
(p + q) : p resonance by a planet 



in^, 

\p + ql flo 



(14) 



where aq and a are the initial and current semi-major axes of 
planet. Supposing a planetesimal gets captured in the 2:1 reso- 
nance with eo - 0, it would have to be dragged to a/oo - 5 to 
reach an eccentricity of 0.9. 

Resonance capture is not guaranteed however. As the res- 
onance of a planet moves in (or out), because of the planet 
migrating, and it encounters a planetesimal, the probability of 
capture depe nds on both eccentricit i es, the migration rate, and 
planet mass (iMustill & Wvattll201(ih . iMustill & WvattI conclude 
resonance capture is guaranteed for low eccentricity and low mi- 
gration rates, impossible at fast migration rates and low eccen- 
tricities, and possible for fast migration rates and high eccen- 
tricities. The transitions between these regimes depend on the 
planet mass. Furthermore, the critical migration rate for c apture 
increases with /', so if the particle survives the passing of the 2:1 
resonance, it may still be captured in another resonance closer to 
the planet with higher j. It is also pointed out that planetesimals 
that are not caught, experience a small jump in e up or down de- 
pending on the migration rate. IMustill & WvattI (1201 Ol) calculate 
capture probabilities as a function of migration rate for differ- 
ent MMRs, migration rates, and planetary masses. For instance, 
a Neptune-mass planet is found to capture all bodies in its 3:2 
resonance for migration rates < 4 x 10"^ AU yr"' . 



Resonance trapping followed by eccentricity-pumping ap- 
pears to be needed to ensure high relative velocities. Another ar- 
gument for an inwardly migrating planet is the fact it can collect 
and carry material inwards. A stationary planet can only scatter 
planetesimals in its direct vicinity. Since the optically thin struc- 
tures are located close to their central stars and have considerable 
mass in small dust grains (see section|2]l this build-up of material 
by an inwardly migrating planet is crucial. 

Assume a planetesimal disk with a mass distribution equal to 
the mass in solids in a minimum mass Solar Nebula (MMSN). 
The surface de nsity of solids in a MMSN disk is given by 
(lHavashilll981h 



[aUj 



-3/2 



gem 



(15) 



So = 7.1 for r < 2.7 AU because ices cannot survive in the hot 
inner regions. Outside the snowline Eo - 30. 

A stationary planet of Neptune-mass on a circular orbit at 5 
AU can sca tter planetesima ls within a region the size of its Hill 
sphere (e.g. lArmitagell2010l) 



Ar - Ru - a 



AW'" 
3M. 



(16) 



The mass of planetesimals within a ring of this size is given by 



Am — InaLpAr. 



(17) 



Using Ep - 30gcm"^, we find Am ^ 4.6M^. 

Assume now that the planet migrated from 30 to 5 AU at 
such a rate, it collected all material in between. We can integrate 
eq. [15] to find Am ^ 46Mq. This is an increase by a factor of 
10 with respect to the mass scattered by a stationary planet at 5 
AU. The difference becomes even bigger when looking at planets 
closer to the star, since Am in eq.[T7]is proportional to a^. 



4. Simulations 

4.1. Introduction 

The scenario for creating the highly-inclined orbits is now as fol- 
lows. A planet with Vesc ^ vk and a planetesimal population start 
out on (nearly) circular orbits and in the orbital plane (figure[3^). 
Due to for example the presence of gas the planet starts migrat- 
ing inwards, and collecting material in mean motion resonances. 
Planetesimals that are stuck in an MMR, have their eccentrici- 
ties raised significantly. Figure |3j) shows the system after some 
migration has taken place. All planetesimal orbits in this figure 
have the same semi-major axis (they are presumed stuck in the 
same MMR) but diff'er in eccentricity. As the planet migrates 
even further, the increasing eccentricities cause the orbits of the 
planetesimals to be planet-crossing (figure|3];). At this point, the 
inclinations are only a couple of degrees. The relative velocity 
between the planet and the planetesimals grows together with the 
eccentricity of the planetesimal orbit and at some point this will 
lead to violent scattering events. As the planet has nearly reached 
the star it will leave behind a scattered population of planetesi- 
mals, which show (depending on the planet's escape velocity) a 
large spread in inclinations, figure |3]l- This population of scat- 
tered planetesimals, over time, will collide and grind down to 
form the small dust that causes the optical depth. In this section 
numerical simulations are presented which show the resonance 
capture, eccentricity-pumping, and scattering of planetesimals 
happening simultaneously. 
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(a) 



pole on 



edge on 



(b) 



(c) 



(d) 





Figure 3. Illustration of the mechanism to reach a planetesimal 
population with high orbits, as explored in this paper (a) both 
the planet orbit (dotted) and planetesimal orbits (solid) start out 
in the orbital plane, (b) As the planet has migrated in a some 
distance, the planetesimals are stuck in an MMR. As they are 
dragged along their orbits become increasingly eccentric, (c) At 
some point the planetesimal orbits become so eccentric, the bod- 
ies become planet-crossing. As the eccentricity rises even fur- 
ther, the relative velocity between the planet and the planetesi- 
mals grows, (d) The high relative velocity compared to the lo- 
cal Keplerian velocity enables the planet to violently scatter the 
planetesimals throughout the system, to orbits with high inclina- 
tions. 



4.2. Setup 

For the integrations of the migrating planet the software pack- 
age Swifter is usecQ. Swifter is written by D.E. Kaufmann 
and based upon the Swift package developed by Levison 
& Duncan. Within Swifter, the user is given the choice 
between 7 integrating schemes, in cluding Wisdom-Holman 
Mapping dWisdom & Holmaiil Il991h . the Regularized Mixed 
Varia ble Symplectic (RMVS) integrator (ILevison & Duncan 
I 994. and the Symp l ectic Massive Body Algo rithm (SyMBA) 
( Duncan etal] 119981: ILevison & DuncanI I200Q) . I n addition 
Swift er also includes a Bulirsch-Stoer integrator (iPress et al.l 
Il992h . In this work SyMBA was used most of the time, since it 
can handle close encounters correctly and is significantly faster 



http://www.boulder.swri.edu/swifter/ 



than the BuUrsch-Stoer method. The results of the SyMBA inte- 
grator are checked by comparing a few cases to a Bulirsch-Stoer 
integration with a very small timestep. 

During the integrations, a planet is pushed to migrate in- 
wards from 30 AU, while a population of 100 massless plan- 
etesimals orbits the star at circular orbits between 10 and 15 AU. 
The planetesimals all start out with e - i - Q, whereas the planet 
is given a small inclination of 0.72°, to ensure <p can take on 
all values between and 27r, and an eccentricity of 0.012. The 
planet mass is varied between 0. 1 and 20 Neptune masses. 

For simplicity, the migration is not done in a self-consistent 
fashion. We want to focus on the effect of the migrating planet 
on the planetesimal population, to see if the resonance cap- 
ture, eccentricity pumping, and scattering can take place. A self- 
consistent migration would vastly increase the complexity of 
the problem. Instead, the planet is pushed within SyMBA to 
have it migrate at a constant rate ^, which is varied between 
7.3 X lO^^AUyr"' and 7.3 x lO^'^AUyr"'. These migration 
rates are based upon esti mates for migration rates and time 
scales for type I migration (Tan aka et al.ll2002l) . type II migra- 
tion (Armi tage 2010), and migration through a planetesimal disk 
dArmitageHSnob . Since the migration in this work is artificial, 
the system does not contain any gas or massive planetesimals 
and, therefore, does not distinguish between the various types of 
migration. Section|6jwill discuss what effect the presence of gas 
or other drivers for migration might have on the results. 

For SyMBA to be able to handle close encounters between 
massless bodies and planets, a Hill-sphere radius, Ru, needs 
to be specified for every planet. This radius determines when 
SyMBA will shorten the integration timestep. Since the Hill- 
sphere radius is proportional to a planet's semi-major axis (eq. 
[TSI l. it will decrease in time for a planet that migrates inwards. 
For simplicity, we will treat Ru as a constant. This results in the 
size of the Hill sphere being overestimated as the planet migrates 
in, slowing down the integrations somewhat. 

The system is evolved long enough for the planet to pass 
through the planetesimal population, and end up close to the 
star, < 1 AU, where it has a negligible effect on the planetesi- 
mals. The timestep for the integrations ranges between 0.5 days 
and 10 ' yr, depending on the smallest perihelion distance in the 
simulation, making sure all orbits are resolved. Within SyMBA, 
the timestep is decreased when a planetesimal enters the Hill 
sphere of the planet. During the simulations, planetesimals are 
removed if they 

- come within a distance r^i„ from the origin. We use rmin = 
0.05 AU. 

- have a > rmax- In our simulations r„inx - 1000 AU is used. 

- come within Rpi. In this work Rp\ - /?Nep is used. 

The choice for r^in comes from the fact that very eccentric or- 
bits with perihelion distances smaller than 0.05 AU are hard to 
resolve without using a very small timestep, and will lead to non- 
physical behaviour Furthermore, bodies that come this close to 
the central star are likely to be evaporated. 



4.3. Resonance Capture 

An example of resonance capture in the simulations is illus- 
trated in fig. |4] The figure shows a Neptune-mass planet be- 
ing pushed artificially to migrate in from 30 AU at a rate of 
7.3x 10"^ AU yr"' . A massless asteroid is orbiting the solar-mass 
star at 13 AU on an initially circular orbit. The asteroid survives 
the passing of the 1:2 (j - 1) MMR but is eventually captured 
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Figure 4. A Neptune-mass planet is migrating in with ^ = 7.3 x 
10"^ AUyr"'. The asteroid survives the passing of the 1:2 reso- 
nance but is captured in the 3:4 resonance and dragged inward. 



by the 3:4 ( j - 3) MMR. For lower migration rates (or higher 
planet masses) the asteroid gets trapped in the 1:2 resonance al- 
ready, while for higher migration rates (or lower planet masses) 
the asteroid only gets captured very close to the planet. 

Figure |5] shows the evolution of a massless body as it is 
dragged along in a migrating planet's 1:2 resonance. The inte- 
gration has been performed using two different integrators, to 
check the consistency of SyMBA for highly eccentric orbits. The 
test particle is captured at f s; 0.13 Myr and eventually lost to the 
central star after about 0.32 Myr. During this time, the lOMNep 
mass planet has travelled ~ 14 AU. The results of the SyMBA 
and Bulirsch-Stoer integrators are very similar, and in both cases 
the massless body is eventually lost to the central star when it 
reaches e x 0.99. The expected increase in eccentricity, as pre- 
dicted by eq. [T4]is also plotted and seen to agree well with the 
simulations. The discrepancy between eq.[T4]and the simulations 
towards higher e comes from the fact eq. [14] is a first-order ap- 
proximation. 

4.4. Inclinations 

Figure shows the evolution of the planetesimal swarm as a 
function of time. The left panels show eccentricity, the right 
panels inclination, and time passes from top to bottom. The 
large circle shows the location of the planet and the small cir- 
cles represent massless planetesimals. Planetesimals above the 
solid lines are on a planet-crossing orbit. The dotted lines show 
the 1:2 and 2:3 MMR, where the 1:2 MMR is furthest from 
the planet. The planet is O.IMncp and migrating at a rate of 
7.3 X 10"^ AUyr^'. The planetesimals start out between 10 and 
15 AU with e - i - 0. After 0.15 Myr (top panels), both reso- 
nances have reached the planetesimal population, and both ap- 
pear unable to capture planetesimals. The passing of the reso- 
nances through the swarm of planetesimals is seen to have a 
very small effect on their orbital parameters. After 0.3 Myr the 
planet itself has migrated through the population of massless 
bodies, and has scattered most of them via direct encounters. 
Since the planetesimal orbits were not excited before the en- 
counter, the relative velocity remained small, and the scattering 
events proved unable to get the planetesimal inclinations higher 
than a couple of degrees. The resulting population, shown after 
0.4 Myr in the bottom panels, is very similar to the one after 0.3 



0.4 



0.2 




50000 100000 150000 200000 250000 300000 350000 
t[yr] 

Figure 5. Evolution of the eccentricity of a massless body stuck 
in the 1:2 internal mean motion resonance of a IOMnbp planet 
migrating inward at a rate of 7 .3 x 10"^ AU yr"' . The timestep for 
both the SyMBA and Bulirsch-Stoer integration was set at 0.5 
days. The dotted line shows the expected increase in eccentricity 
as predicted by equation [14] The test particle is captured at f s; 
0.13 Myr and eventually lost to the central star after about 0.32 
Myr. 



Myr, i.e., the planet has little effect on the bodies once it has mi- 
grated inside the population. During this simulation, 2% of the 
planetesimals is accreted by the migrating planet. 

If we increase the planet mass, this will lead to a high cap- 
ture probability. Figure |6j) shows a similar calculation, but for 
a planet mass of 2MNep- After 0.15 Myr the effect of the en- 
hanced planet mass are already visible as the passing of the 1:2 
resonance raises the eccentricities of the planetesimals slightly, 
while the 2:3 resonance appears to be able to capture bodies, and 
subsequently drag them along. This becomes even more visible 
after 0.3 Myr, where a significant fraction of the planetesimals 
are stuck in the 2:3 MMR. As these bodies are dragged along, 
their eccentricities are raised as described by eq. [14] As the ec- 
centricity of a planetesimal goes up, it will eventually reach a 
planet-crossing orbit, allowing it to have a close encounter with 
the planet. This has already happened for the bodies no longer 
in the 2:3 resonance. Since the close encounters happen when 
the planetesimals are on eccentric orbits, the relative velocities 
involved are higher, and the inclinations resulting from the scat- 
tering events are much greater than for figure During this 
simulation, 1 % of the bodies were accreted by the central star, 
and 1 % by the migrating planet. 

If we increase the planet mass even further, the capture prob- 
ability reaches almost unity, and the vast majority of planetesi- 
mals get stuck in a resonance far away from the planet. As they 
are dragged along and their eccentricities go up, they are pro- 
tected from a close encounter with the planet for some time, 
since they are so far away. It turns out, almost all these bodies 
will reach e » 0.9, and are accreted by the central star before 
they have time to be scattered by the migrating planet. An ex- 
ample of this scenario is shown in figure [5] This means, that in 
the systems where you would expect to find the highest relative 
velocities (and therefore the highest inclinations) the majority of 
the planetesimals are actually lost to the central object. A way to 
get around this, might be to put an additional, stationary planet 
close to the central star, which will scatter the highly excited 
planetesimals before they come too close to the central object. 
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Figure 6. Evolution of the eccentricities (left) and inclinations (right) in a planetesimal swarm as a planet (large circle) migrates 
inwards from 30 AU at a rate of 7.3 x 10"^ AUyr"'. The planet mass used are O.IMnbp for (a), and 2MNep for (b). Planetesimals 
(small circles) above the solid lines are on a planet crossing orbit, and the dotted lines show the location of the 1:2 and 2:3 MMR. 



This planet could not be too massive, as this would lead to it 
ejecting the planetesimals on hyperbolic orbits. In this study we 
have restricted ourselves to single-planet systems, but the study 
of stellar systems with multiple planets, one of which is migrat- 
ing, might be very interesting indeed. Another way to protect the 
planetesimals from accretion by the central star would be to in- 
troduce a gas drag, which would dampen the eccentricities. The 
effect of gas of the results is further discussed in section |6] 



From section |43] we know that the capture probability also 
depends on the migration rate of the planet. We have therefore 
conducted a small parameter study, where the migration rate and 
planet mass have been varied. In total we carried out 12 sim- 
ulations with the planet-mass ranging from 0.1 - 20MNep and 
^ between 7.3 x lO""* - 7.3 x lO-*" AUyr"'. Table □ shows the 
characteristics of planetesimal populations after migration has 
taken place, as well as the fraction of planetesimals which have 
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Table 1. Orbital parameters of the resulting planetesimal population for a number of simulations. 



Sim 


^[AUyr"'] 


Mpl/MNep 


A'ej[%] 


A'acc[%] 


<fl>[AU] 


{e) 


<0[°] 


^max [ ] 


SO.l 


7.3 X 10"^ 


0.1 








9.32 


0.255 


7.39 


16.16 


SI 


7.3 X 10"* 


1 





100 


- 


- 


- 


- 


MO.l 


7.3 X 10"' 


0.1 








12.72 


0.047 


1.39 


3.98 


Ml 


7.3 X 10"' 


1 








10.51 


0.311 


8.13 


19.00 


M2 


7.3 X 10"' 


2 





1 


8.61 


0.538 


14.61 


35.25 


MIO 


7.3 X 10"' 


10 





100 










M20 


7.3 X 10"' 


20 





100 










FO.l 


7.3 X 10"" 


0.1 








12.61 


0.016 


0.29 


1.84 


Fl 


7.3 X 10"" 


1 








12.59 


0.082 


1.84 


5.92 


F2 


7.3 X 10"" 


2 








13.23 


0.113 


2.50 


9.10 


FIO 


7.3 X 10"" 


10 








12.30 


0.442 


9.60 


27.86 


F20 


7.3 X 10"" 


20 








18.79 


0.408 


9.42 


47.23 



Notes. A'ej and W^cc give the percentage of planetesimals that were lost during the integration resulting from respectively ejection from the system, or 
accretion by the central star. Inclinations are given in degrees. Simulation names consist of one letter describing the migration (F=Fast, M=Medium, 
S=Slaw) and a number representing the planet mass. 



been lost during the integration. Both higher planet mass and 
lower migration rate appear to ensure high average inclinations. 
This is the case because; 1) massive planets are more capable of 
scattering planetesimals into highly-inclined orbits (see section 
13.2b . and 2) relatively massive planets and low migration rates 
increase the probability of resonance capture, allowing for more 
eccentricity-pumping and a higher relative velocity between the 
planetesimal and the planet during the close encounter. However, 
if the combination of migration rate and planet mass causes the 
planetesimals to get stuck in a resonance too far, they are likely 
to be lost to the central star before having a close encounter with 
the migrating planet. 



5. Dust Production 

Section|4]has left us with populations of planetesimals and con- 
straints on their orbital parameters. Section 15. II will calculate 
the time scales on which these planetesimals should collide and 
form dust, and section |S!2] will estimate the total mass and sur- 
face area of the dust size distribution resulting from the colh- 
sional cascade. 



5. 1 . Collision Times 

Suppose we have planetesimals of radius R, mass m and a 
geometrical cross-section for collisions of ctcoh = A-nR^, flying 
around in some volume V. We can define the sweeping time as 



V 



VconC coll 



(18) 



where Vcon is the collision velocity between the planetesimals. 
Since the populations from table [T] show pretty eccentric orbits, 
the collision velocity will vary a lot between colhsions, but for 
the purpose of this calculation we estimate it as 



Vcoll = rjVK «fl» , 



(19) 



with T] between 0.5 and 2. For two nearly circular orbits rj is 
determined by the vertica l component of the individual velocities 
dPominik & Decinl2003b . But since the orbits in our populations 
are eccentric and randomly oriented, 77 will lie closer to V2. 



If we assume the volume is wedge-shaped and lies between 
an inner and an outer radius, rin and rom, and has a normalized 
height of li = H/r, we can write it as 



An 



-hiri 



(20) 



where h - tan/^ax- 

We have to think carefully about and rout- One might in- 
sert the minimum and maximum values of the semi-major axis 
found in the resulting populations, but this is not correct since 
most comets are on pretty eccentric orbits. On the other hand, 
choosing rout to be equal to flmax(l + ^max) will increase the vol- 
ume drastically even though only one planetesimal is able to get 
to rout. It makes more sense to define something like 

rou, =«fl> + fr„)(l+(e)), (21) 
rin = ({a) - CT„)(1 - (e)), (22) 

where {a} and (e) can be taken from table [T] We approximate 
~ V(fl). which is found to work well for most populations in 
tables 

The collision time for the comets will become 



Wi - ^, 



(23) 



with the number of bodies being constrained by the total mass 
of the population M 



m \ R. 



M 



(24) 



where p denotes the density of the planetesimal material. We can 
now calculate the collision time scale for diff'erent populations. 
A promising case for creating a halo-like structure appears to 
be M2, which has high inclinations but is still located relatively 
close to the star. From table[T]we find (a) - 8.61 AU, <e) = 0.538 
and /max - 35.25°. For this specific configuration of comets, the 
collision time becomes 



(M2) _ 3 7 
'coll -J-'^ 



1 1.41/ \km/ 



P 



M 



0.5pe/\50Mg 



yr. (25) 



The collision time scale for the population in M2 is compa- 
rable to the migration time itself. This means a number of as- 
teroids might suffer collisions during the migration itself, which 
we have not considered. The collision time is also an indication 
of how long the dust can survive in the system, as new dust is 
continuously formed during this period. 
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5.2. Dust Size Distribution 

When planetesimals collide at sufficiently high velocities, they 
will fragment and create smaller bodies. The distributions of 
these bodies has the shape of a power law 



f(m) oc m 



(26) 



with q - 11/6 (lDohnanvilll968h . Assuming spherical particles 
we can rewrite this as a size distribution 



f(R) = M^, 



(27) 



withy = 2-3q = -3.5. The proportionality factor /r depends on 
the material properties of the initial population and is calculated 
bv lDomin ik & Decin (2003) as 



2m^^ 
\(4;r/3)peo V eo 



with 



and 



y + 1 



''coll 

4S 



"coll 



"coll 



1652 2S 



- 1, 



(28) 



(29) 



(30) 



where S is the binding energy of the material. F or asteroid type 
bodies 5 = 10^ erg g ' dPominik & Decinll2003h . 



5.3. Total Mass and Optical Depth 

The size distribution can be integrated to yield the total mass 



R^pf(R)dR, 



(31) 



in dust grains with radii between R^i^ and R^ilx- 

The smallest size, Rmin, is determined by comparing the 
gravitational force acting on the dust grain with the radiation 
pressure 



3Zy* Qpr 

l6ncGM,Rp' 



(32) 



with c the speed of light and Qp^ the radiation pressure coefficient 

Gpr = gabs + Qsc-Jl - (COS a}), (33) 

where a denotes the scattering angle. F or isotropic forward scat- 
tering, (cos a) - I dBums et al.lll979l) . Particles with /3 > 1/2 
are ejected from the system. 

Assuming 2pr = 1, a dust grain with p - 2.76 gem"-' or- 
biting a sun-like star reaches the critical value of/? = 1/2 for a 
radius of about 0.4 /zm. From here on on we will assume lOfim 
to be the maximum size of a dust particle still contributing sig- 
nificantly to the optical depth. 

Again considering M2, we have Vcoii - 14.4 km s"^ Together 
with S - 10^ erg g"' this gives eo = 1.90 X 10^ If we take the 
total planetesimal mass to be 50M^, and the material density 
of the bodies to be 2.76 g cm"^, we find the total mass in dust 
grains between 0.4 and lO/vm to be Mdust = 8.7 x IO^'^Mq. 
This mass compares well to the mass in small dust grains quoted 
in section|2] but depends heavily on parameters such as the total 



planetesimal mass, and is only intended as an order of magnitude 
estimate. 

The optical depth along a line of sight dr can be written as 



dr — -Kpdr — ncrdz. 



(34) 



with K the opacity and n the number density of the particles re- 
sponsible for the absorption. Assuming the dust is evenly dis- 
tributed we can rewrite this as 



, cr (Just , 
dr — ^ dr. 



(35) 



where V is the volume in which the dust is located, and o-^ust is 
the total combined cross-section for radiation of the dust grains 



O" dust 



r 



nR^f(R)dR. 



(36) 



Assuming the dust traces the planetesimal population, i.e., de- 
structive collisions happen throughout V, we can integrate eq. 
|35]find 



2.76 g cm 



-3 



M 

50Mffl 



(37) 



The radial optical depth is indeed < 1 . Again, the derived opti- 
cal depth depends a lot on the initial parameters. Since we have 
assumed the dust density is constant in the volume V, the radial 
optical depth is almost independent of the angle away from the 
midplane. 

The optical depth could be increased by making the planetes- 
imal population more massive, or changing their size. This will 
however also increase Mjust- 

A better way to increase r might be to confine the volume in 
which the dust is located. This can be done in two ways; either 
moving the entire planetesimal population closer to the star, or 
by realizing the planetesimals have pretty eccentric orbits and 
mo st destructiv e collisions might therefore happen near periap- 
sis. lWvattetal.l(l2010) showed, using a Monte-Carlo simulation, 
that a population of bodies with identical semi-major axes and 
eccentricities, but random mean longitudes, arguments of peri- 
centre and longitudes of ascending nodes, 90% of destructive 
collisions happen within a radius rja - (1 - e^)l(\ - 0.72e), 
which might be a hint in this direction, even though the bodies 
in the populations in table[T]are distributed both in a and in e. 

6. Discussion 

Comparing the SEDs of several YSOs to radiative transfer mod- 
els and near-IR interferometric visibilities appears to indicate the 
presence of an optically thin dust structure with a large scale 
height (section |2]i. We have argued that an inwardly migrat- 
ing planet could explain the observed structure. The scenario is 
that an inwardly migrating planet collects planetesimals in mean 
motion resonances, raises their eccentricities, and subsequently 
scatters them to orbits with greater inclinations (fig. |3). Over 
time, the planetesimals collide and form small dust in a colli- 
sional cascade. 



6.1. The role of eccentricity pumping 

An important result of this study is that scattering planetesimals 
efficiently into highly-inclined orbits requires the relative veloc- 
ity between the planet and the planetesimal to be high, of order 
of the local Keplerian speed. This is the reason why in our setup 
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the sequence of eccentricity pumping in resonances followed by 
dynamical scattering is the most successful strategy. 

The resonance capture probability of a migrating planet is 
a function of migration rate, planet mass, an d distance to the 
star, a nd increases for more massive planets (iMustill & WvattI 
I2OIOI) . This puts constraints on the planet mass and migration 
rate. The process of eccentricity-pumping is well-defined by 
iMalhotra (1995); Wyatt (2003) and eq. M\ 

can be used to es- 
timate the distance over which a particle has to be dragged in 
a particular MMR for its orbit to reach a certain eccentricity. 
However, if a planetesimal gets stuck in a resonance that is too 
far from the planet, it is likely to be accreted by the central star 
after its eccentricity approaches unity. This fate could possibly 
be averted by introducing gas drag, or a second, stationary planet 
close to the central star. In this case, the migrating planet would 
excite and push the planetesimals inwards, while the stationary 
planet would cause the scattering events. The mass of this second 
planet should not be too high, as this would lead to the ejection 
of the excited planetesimals. Systems with Neptune-mass plan- 
ets within an AU of their parent stars are not uncommon (e.g. 
[Bouchv et al. 2009), and adding such a planet to the simulations 
presented in this work could well lead to even higher planetesi- 
mal inclinationS;_^ 

IYu & Tremaind (1200 ih study resonant capture by a Jupiter- 
mass planet migrating inward using analytic arguments and 
three-body integrations. They find that in typical systems, the 
resonance capture and subsequent eccentricity-pumping will 
most likely result in accretion of the test particle by the cen- 
tral star. This fate can be avoided by having a close encounter 
with the migrating planet, which eventually leads to an escape 
from the system. In the study presented here, systems with very 
massive planets show very similar behaviour. 

6.2. Migration 

Our numerical simulations show that the basic scenario does 
work well for a certain range in planet mass and migration 
rate. The migration rates used are consistent with migration time 
scales for type 1, type II, and migration through a planetesimal 
disk. However, the migration itself is not done self-consistently 
in our simulations. Because of this we have not taken into ac- 
count the effect of the driver for the migration on the planetesi- 
mal population. 

One possible driver is gas. A planet in a gaseous disk can ex- 
change energy and angular momentum with the gas. If the planet 
is unable to significantly alter the disk structure it is referred to as 
type I. In this scenario the exchange of angular momentum can 
be written as the sum of the torques exerted by the Lindblad res- 
onances of the planet on the disk (Goldreich & Tremaine.,1979^ 
Il980l) . The direction of the migration depends on the size of the 
interior and exterior torques and is found to be inward for an 
isothermal disk (Ward 1997). Paardekooper & Mellema (200^ 
simulated non-isothermal disks and showed the direction of mi- 
gration dep ends on the ability o f the disk to radiate away gen- 
erated heat. lTanaka et al.l (l2002h calculated the sizes the torques 
for isothermal disks and derived migration time scales of the or- 
der of 1 Myr. The migration time scale is inversely proportional 
to planet mass, massive planets migrate faster (I Armitageil20 1 Ol) . 

A slightly more massive planet will be able to open up a gap 
in the disk. In this case the co-rotating torque disappears. The 
migration rate for type II migration is equal to the gas-inflow 
velo city and found to be independent of planet mass ( Armitage 
I2OIO ). Type II migration is generally faster than type I and at 5 
AU a typical disk results in a migration time scale of 10^ yr. 



Gas is not necessarily needed for migration though. A planet 
embedded in a planetesimal belt can exchange angular momen- 
tum with these bodies during scattering events. If the planet scat- 
ters bodies inward and outward at an equal rate, the net change 
in the angular momentum of the planet will be zero. If however, 
the planet usually scatters bodies to larger orbits, it w ill slowly 
loose angular momentum and migrate towards the star (llda et al.l 
2000). For the change in angular momentum to have a signifi- 
cant effect on the planet's orbital radius, the mass of the planet 
has to be lower or equal to the total mass of scattered bodies. 
Also, for the migration not to ha lt the pla n et ha s to encounter 
fresh planetesimals continuously, llda et al. I (120001) calculate the 
migration rate to be 



da 
dt 



pi 



(38) 



with the disk mass equal to Mdisk - ^^Opi- iKrrsh et al] (l2009l) 
study the orbital evolution of an isolated planet in a planetesimal 
belt and confirm this result, adding that for typical planetesimal 
surface density profiles the direction of migration will generally 
be inward. 

It is stressed that the planetesimals that are having their ec- 
centricities raised in an interior resonance are not the ones re- 
sponsible for the migration. Only when these bodies reach very 
eccentric orbits and suffer close encounters with the planet are 
they able to contribute to the planetary migration. 

Aside from planet-disk interactions, the semi-major axis of 
a planet can also be chang ed significantly as a re sult of inter- 
actions with other planets. [Raymond et al.l (1201 Ih numerically 
investigated systems with an initially unstable configuration of 
giant planets within a debris disk. It was found that the dynami- 
cal instability could trigger a coUisional cascade. 

6.3. The influence of gas in fhe disk 

While planetesimal scattering alone can make a planet migrate, 
the preferred scenario may be to have the planet migrate in the 
phase when the disk is still gas-rich, while our simulations with 
forced migration have ignored the effects of gas on the planetes- 
imals. If we attribute the planet's radial motion to type I, type 
II, or planetesimal-driven migration in a gas-rich disk, we need 
to consider the effects of the gas during the migration and scat- 
tering phases. The presence of gas will cause aerodynamic gas 
drag. This drag will have two main implications; it will dampen 
the eccentricity and inclination of a planetesimal, and, it will 
have an effect of the probability of capturing a planetesimal in 
an MMR. Since the capturing of planetesimals followed by the 
eccentricity-pumping plays a vital role in this work, the effects 
of gas on these two mechanisms have to be addressed. 

The aerodynam ic gas drag time scale is given by 
lAdachietal.l(fl976l) 



8pR 



3Cz3PgasVK 



(39) 



with pgas and p the density of the gas and planetesimal mate- 
rial respectively. The drag coefficient Co is of order unity, and 
a non-linear function of the planetesimal's radius R and relative 
velocity to the gas. From the equation it is clear that gas drag will 
predominantly effect smaller bodies at relatively small distances 
to the star. 
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ICapobianco et alJ (1201 ih study planetesimal driven migra- 
tion in a gas disk and find the eccentricity of a planetesimal is 
dampened by a factor of order itself on a time scale of 



2.3 X 10V/l-0\/ R 



f(a, e) \ 



l.UW R \l a 
IT/lkm/lAU 



13/4 



(40) 



where = 1 corresponds to a MMSN-disk and f{a, e) is given 
by 



2/3 



X 



(41) 



with;^' = Rn/api and 770 - 1.95x10"-' AU"'^^. The damping time 
scale is smallest for small bodies since they couple to the gas bet- 
ter. For a larger planetesimal of 10 km with a semi-major axis of 
10 AU, the damping time scale i s ~ 10^~^ yr and the system ap- 
proaches the gas-free case (Cap obianco et al.ll201 ll) . It appears 
that the eccentricity-damping plays a small role for planetesi- 
mals of sizes > 1 - 10 km. For smaller planetesimals or higher 
gas densities, gas drag will play an important role. From eqs.l40l 
and l4n it is obvious gas drag will affect eccentric orbits more 
than circular ones. In such a regime, eccentricity-damping could 
therefore protect planetesimals against accretion by the central 
star. This would allow for them to be scattered by the migrating 
planet, even if they are initially captured in a resonance far from 
the planet. 

The coupling of planetesimals to the gas has an effect on 
the p robability of resonance capture as well. ICapobianco et alj 
(1201 1) define a radius /?trap for which half of the planetesimals 
are captured. For a MMSN disk this radius corresponds to 



— (g)(^r(sr(^r""--> 

with Cd the drag coefficient of order unity. For a Neptune-mass 
planet at 15 AU, ^drag ~ 0.01 km. Assuming the planetesimals 
are of km-size, the effect of gas on the capture probability can 
thus be neglected. 

TFogg & Nelsoiil (l2005h have studied the effect of a giant 
planet migrating in on the process of terrestrial planet formation, 
using N-body simulations. The inward migration is implemented 
artificially with a prescription that is chosen to match type II mi- 
gration. The disk is assumed to be gas-rich. Consequently, the 
gas drag causes radial drift of planetesimals, an d a damping of 
eccentricities and inclinations. iFogg & Nelsori find that the in- 
ward migration causes the planet to capture small bodies in its 
inner resonances, after which the orbits of these bodies are ex- 
cited. The excitation of the orbits, followed by direct scatter- 
ing by the planet, causes the formation of an exterior scattered 
disk of planetesimals, which looks very similar to the resulting 
populations in figure|6] Depending on the initial disk properties, 
the exterior disk is found to contain 30-70% of the inner disk 
material. The gas damping is found to play an important role 
in the inner reg ions of the disk. For increasingly mature disks, 
iFogg & NelsonI find a wider spread in the orbital parameters of 
the scattered disk, caused by the increasingly smaller effect of 
gas damping. 

In a somewhat similar study, with the focus of fo rming Earth- 
like pl anets in the habitable zone, is conducted bV iMandell et al.l 
(l2007h . They find si milar results for an in wardly migrating 
Jupiter-mass planet as Fogg & N elsonI (l2005l) . including a scat- 
tered exterior disk. They also show that a second, stationary gi- 
ant planet outside of the migrating one, will remove this scat- 
tered disk by ejecting the bodies from the system. A smaller or 



vanishing gas drag leads to more material ending up in the outer 
disk. The spread in orbital parameters of the bodies in the scat- 
tered disk also increases when less gas is present. 

The extended scattered disks in these bodies show similar 
characteristics as the populations found in table [Tj when gas 
does not play a big role. However, these studies focus on gi- 
ant planet migration, in relatively young disks, thus focussing 
on low-velocity collisions between the planetesimals that allow 
them to grow and form planet embryos and Earth-like planets. 
In the study presented in this paper, we are instead focussing on 
less massive planets, and interested in the orbital parameters of 
planetesimals in the scattered disk. 

After the migration has taken place, the planetesimals collide 
on a certain time scale to form dust. The dust distribution can be 
calculated under some approximations and the total mass and 
radial optical depth are found to be consistent with the observed 
haloes for certain values of variable parameters such as the to- 
tal mass in the disk, the initial size of the planetesimals, and the 
planetesimal material density. The coUisional model used in sec- 
tion |5] and the estimates of the optical depth and total mass are 
very simple and have their limitations. For instance, we assume 
the dust to trace the comet population, use simple geometrical 
estimates to define the volume the planetesimals occupy, and ne- 
glect the effect of any gas that might be present. Because of the 
nature of the discussed systems, including the spread in semi- 
major axes, eccentricities and inclinations, it is difficult to avoid 
these simplifications without having to make a detailed coUi- 
sional and radiative transfer model, which is beyond the scope 
of this paper. If gas is present in the inner regions of the system, 
it will have a significant effect on the small dust grains since they 
couple to the gas on short time scales. Gas which is confined to 
the disk might rapidly remove small dust grains from highly- 
inclined orbits and reduce the scale height of the optically thin 
structure. Since we are dealing with pre-transitional disks, it is 
very well possible the inner regions of the system do not contain 
any gas at all. 

7. Conclusions 

The goal of this work was to see if an inwardly migrating planet 
could be the cause of the observed optically thin halo-like struc- 
tures. The main findings are: 

- An inwardly migrating planet is capable of capturing and 
transporting material inward in mean motion resonances. 
The probability of capture depends on planet mass and mi- 
gration rate. Larger planet masses and lower migration rates 
increase the capture probability. The presence of gas on the 
capture probability can be ignored for km-sized planetesi- 
mals. 

- When a trapped planetesimal is pushed inwards, its ec- 
centricity will rise rapidly and in agreement with equation 
[T4I During this phase the inclination does not reach values 
greater than a couple of degrees. The growing eccentricity 
will eventually put the planetesimal on a planet-crossing or- 
bit while increasing the relative velocity between the two. 

- Simulations show that the inward migration of a planet can 
lead to a population of planetesimals on highly-inclined 
orbits. The average inclination of the remaining popula- 
tion increases for higher planet masses and lower migration 
rates. Resonance capture and subsequent eccentricity pump- 
ing plays an important role in this outcome. 

- The high planetesimal inclinations in these systems result 
from direct scattering events, and reach significantly higher 



13 



S. Krijt and C. Dominik: A dynamical model for halo-like dust clouds in transitional disks 



values if the planetesimal orbits are excited prior to this 
event. 

- In single-planet systems where planetesimals get stuck in an 
MMR far from the migrating planet, the vast majority of 
planetesimals is lost to the central object. This fate could 
perhaps be averted by introducing a second, stationary planet 
close to the central star. This planet should not be very mas- 
sive, as close encounters would then lead to ejection of the 
excited planetesimals. 

- The resulting population of planetesimals on highly-inclined 
orbits have coUisional time scales of the order of 10^"^ yr 
depending on the properties of the individual planetesimals 
and their combined mass. Since the relative velocities be- 
tween the planetesimals are high (~ 10 km s"') the collisions 
will be destructive and will start a coUisional cascade that 
will grind the bodies down to dust grains. The total mass and 
optical depth of this dust can be estimated and found to be of 
the same order as the total mass and optical depth observed 
in various systems. 

Future studies should focus on self-consistent migration and 
a more careful description of the colUsional cascade and opti- 
cal depth determination to see if this mechanism can work and 
constrain the parameter-space further. With these future studies, 
it might be possible to use the optically thin structures as sign- 
posts for planet formation and migration. The presence of an 
optically thin halo-Uke structure could be induced from an SED, 
without having to resolve the system. The properties of the halo 
could then put constraints on the planet mass, migration rate, and 
distance it has travelled. 

Acknowledgements. The authors would Uke to thank M. Wyatt for useful discus- 
sions, and the referee, P. J. Armitage, whose comments have led to a significant 
improvement of the manuscript. 

References 

Adachi, 1., Hayashi, C, & Nakazawa, K. 1976, Progress of Theoretical Physics, 
56, 1756 

Armitage, P. J. 20 1 0, Astrophysics of Planet Formation, ed. Armitage, P. J. 

Benisty, M., Natta, A., Lsella, A., et al. 2010, A&A, 511, A74+ 

Bouchy, F., Mayor, M., Lovis, C, et al. 2009, A&A, 496, 527 

Brown, M. E. 2001, AJ, 121, 2804 

Burns, J. A., Lamy, P L., & Soter, S. 1979, Icarus, 40, 1 

Capobianco, C. C, Duncan, M., & Levison, H. F 2011, Icarus, 211, 819 

Chiang, E. 1. & Goldreich, R 1997, ApJ, 490, 368 

Dohnanyi, J. S. 1968, in lAU Symposium, Vol. 33, Physics and Dynamics of 

Meteors, ed. L. Kresak & P. M. Millman, 486—1- 
Dominik, C. & Decin, G. 2003, ApJ, 598, 626 
Dullemond, C. R, Dominik, C, & Natta, A. 2001, ApJ, 560, 957 
Dullemond, C. P, Hollenbach, D., Kamp, I., & D'Alessio, P. 2007, Protostars 

and Planets V, 555 
Duncan, M. J., Levison, H. F, & Lee, M. H. 1998, AJ, 1 16, 2067 
Espaillat, C, Calvet, N., D'Alessio, R, et al. 2007, ApJ, 670, L135 
Espaillat, C, D'Alessio, R, Hemandez, J., et al. 2010, ApJ, 717, 441 
Fogg, M. J. & Nelson, R. R 2005, A&A, 441, 791 
Goldreich, R & Tremaine, S. 1979, ApJ, 233, 857 
Goldreich, R & Tremaine, S. 1980, ApJ, 241, 425 
Gulbis, A. A. S., Elliot, J. L., Adams, E. R., et al. 2010, AJ, 140, 350 
Hahn, J. M. & Malhotra, R. 2005, AJ, 130, 2392 
Haisch, Jr, K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 
Hayashi, C. 1981, Progress of Theoretical Physics Supplement, 70, 35 
Hillenbrand, L. A., Strom, S. E., Vrba, F. J., & Keene, J. 1992, ApJ, 397, 613 
Houk, N. 1978, Michigan catalogue of two-dimensional spectral types for the 

HD stars, ed. Houk, N. 
Ida, S., Bryden, G., Lin, D. N. C, & Tanaka, H. 2000, ApJ, 534, 428 
Jang-Condell, H. & Kuchner, M. J. 2010, ApJ, 714, L142 
Jewitt, D. C. & Luu, J. X. 2000, Protostars and Planets IV, 1201 
Kenyon, S. J. & Hartmann, L. 1987, ApJ, 323, 714 

Kirsh, D. R., Duncan, M., Brasser, R., & Levison, H. F. 2009, Icarus, 199, 197 
Levison, H. F. & Duncan, M. J. 1994, Icarus, 108, 18 



Levison, H. E & Duncan, M. J. 2000, AJ, 120, 2117 

Li, J., Zhou, L., & Sun, Y. 2008, Chinese Astron. Astrophys., 32, 409 

Malhotra, R. 1993, Nature, 365, 819 

Malhotra, R. 1995, AJ, 110, 420 

Malhotra, R., Duncan, M. J., & Levison, H. F. 2000, Protostars and Planets IV, 

1231 

Mandell, A. M., Raymond, S. N., & Sigurdsson, S. 2007, ApJ, 660, 823 

Min, M., Dullemond, C. R, Dominik, C, de Koter, A., & Hovenier, J. W. 2009, 

A&A, 497, 155 
Morbidelli, A. & Valsecchi, G. B. 1997, Icarus, 128, 464 
Mulders, G. D., Dominik, C, & Min, M. 2010, A&A, 512, Al 1 + 
Murray, C. D. & Dermott, S. F. 1999, Solar system dynamics, ed. Murray, 

C. D. & Dermott, S. F 
Mustill, A. J. & Wyatt, M. C. 2010, ArXiv e-prints 
Paardekooper, S. & Mellema, G. 2006, A&A, 459, L17 

Press, W. H., Teukolsky, S. A., Vetterling, W. T, & Flannery, B. R 1992, 
Numerical recipes in FORTRAN. The art of scientific computing, ed. Press, 
W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 
Raymond, S. N., Armitage, P. J., Moro-Martfn, A., et al. 2011, ArXiv e-prints 
Safronov, V. S. 1966, Soviet Ast., 9, 987 
Tanaka, H., Takeuchi, T, & Ward, W. R. 2002, ApJ, 565, 1257 
Thalmann, C, Grady, C. A., Goto, M., et al. 2010, ApJ, 718, L87 
Verhoefl', A. P., Min, M., Panfin, E., et al. 201 1, ArXiv e-prints 
Waelkens, C, Waters, L. B. E M., de Graauw, M. S., et al. 1996, A&A, 315, 
L245 

Ward, W. R. 1997, Icarus, 126, 261 
Weidenschilling, S. J. 1975, AJ, 80, 145 

Wilner, D. J., Holman, M. J., Kuchner, M. J., & Ho, P T. P 2002, ApJ, 569, 
L115 

Wisdom, J. & Holman, M. 1991, AJ, 102, 1528 
Wyatt, M. C. 2003, ApJ, 598, 1321 
Wyatt, M. C. 2008, ARA&A, 46, 339 

Wyatt, M. C, Booth, M., Payne, M. J., & Churcher, L. J. 2010, MNRAS, 402, 
657 

Yu, Q. & Tremaine, S. 2001, AJ, 121, 1736 



14 



